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Abstract 

The mass matrix forms of quarks and leptons are discussed in theory with permutation 
flavor symmetry. The structure of scalar potential is analyzed in case that electroweak 
doublet Higgs fields have non-trivial flavor symmetry charges. We find that realistic forms 
of mass matrices are obtained dynamically in the vacuum of the theory, where some 
of Higgs bosons have vanishing expectation values which lead to vanishing elements in 
quark and lepton mass matrices. Mass textures are realized in the true vacuum and their 
positions are controlled by flavor symmetry. An interesting point is that, due to the flavor 
group structure, the up and down quark mass matrices are automatically made different 
in the vacuum, which lead to non-vanishing generation mixing. It is also discussed that 
flavor symmetry is needed to be broken in order not to have too light scalars. The lower 
bounds of Higgs masses are derived from the experimental data of flavor-changing rare 
processes such as the neutral K meson mixing. 



1 Introduction 



Flavor symmetry is expected to be a clue to understand the masses and mixing angles of 
quarks and leptons. It reduces the number of free parameters in Yukawa couplings, and some 
testable predictions of masses and mixing angles generally follow. The discrete non-Abelian 
symmetries contain the symmetry groups Sn, the dihedral groups Dn, the binary dihedral 
(quaternion) groups Qn and the tetrahedral groups A^. Some predictive models with discrete 
flavor symmetries have been explored by many authors [l]-[8]. 

The simplest discrete non-Abelian group is S^. The Higgs sector has a rich structure in 
theory with the 6*3 flavor symmetry, e.g. there are two possibilities that SU{2) doublet Higgs 
fields have trivial and non-trivial flavor symmetry charges. The former case has been discussed 
intensively in the literature (see for example Ref. [2|). The latter case is also expected to have 
characteristic phenomenology where Yukawa couplings of quarks and leptons are described by 
renormalizable operators. 

In this paper, we discuss the mass matrix forms of quarks and leptons in the case that 
SU (2) doublet Higgs fields have non-trivial S3 flavor charges. We perform the analysis of Higgs 
potential at the electroweak scale and examine whether some of Higgs bosons have vanishing 
expectation values. In this case, vanishing elements (texture zeros) of fermion mass matrices 
are obtained dynamically in the vacuum of the theory, and their positions are controlled by 
flavor symmetry. Such zeros in the mass matrices have often been assumed by hand [9]-[17j. 
unlike our scheme. An interesting point of the scheme is that, due to the S3 group structure, 
the up and down quark mass matrices automatically have different forms at the vacuum, which 
lead to non-vanishing generation mixing. The suppression of flavor-changing neutral currents 
(FCNC) mediated by multiple Higgs fields will also be discussed. 

This paper is organized as follows. In Section 2, we present some fundamentals of the 5*3 
group and the 5*3 symmetry invariant form of Yukawa couplings (mass matrices) in the super- 
symmetric case. In Sections 3 and 4, the symmetry-invariant Higgs potential is constructed 
and its structure is analyzed, especially focusing on the conditions for vanishing vacuum expec- 
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1a 1a Is 2 
2 2 2 2 + 1a + 1s 

Table 1: 5*3 tensor products. 



tation values (VEVs). In Section 5, we show that fermion mass matrices favored by the recent 
experimental data are highly limited in our framework. In Sections 6 and 7, several other 
phenomenological issues, i.e. the Higgs spectrum with flavor symmetry breaking effects and 
tree-level FCNC processes, are investigated. Section 8 is devoted to summarizing our results. 

2 5*3 invariant mass matrix 

In this section, the S'3-invariant forms of mass matrices are presented. For more details of the 
5*3 group structure and symmetry-invariant forms, see, for example, Ref. |2]. The S3 group has 
three irreducible representations; two singlets and one doublet, which we denote throughout 
this paper as Ig (singlet), 1a (pseudo singlet) and 2 (doublet), respectively. The decomposition 
of tensor products is shown in Table [H A non-trivial product is 2x2 = 2 + 1a + Is- For two 
doublets ip = {ipi,ip2)'^ and = (0i,02)^, that is explicitly given by 

= (^f^'^ +(^t^^_^t^,),^ + (^t^^ + ^t^2)is , (2.1) 

where the suffixes in the right hand side of the equation denote 6*3 representations. Another 
form of product is written by using ipc = ("^25^1)"^ which behaves as a doublet. The tensor 
product of ipc with another doublet (j) becomes 

^CX0 = (^Jlj^J^ +(V^102-^20l)u + (^102 + ^2</>l)ls • (2.2) 

Since this latter form of product does not contain any complex conjugates, it is useful for de- 
scribing, e.g. Majorana masses for neutrinos and holomorphic terms in supersymmetric theory. 
It is mentioned that the two types of S3 invariants given above correspond to two invariant 
tensors of SU{2), which contains 5*3 as a subgroup. 



Let us discuss quark/lepton mass matrices, where both left-handed and right-handed fermions 
transform under a single symmetry. We suppose that two of three generations belong to 
5*3 doublets and the others are singlets. In this paper, we consider supersymmetric theory and 
then matter superfields of first two generations ^I/li 2 ^^1,2 treated as S3 doublets, into 
which left and right-handed fermions are embedded as ip^^ ^ and ^ , respectively. As for 
SU (2) weak doublet Higgses, all three types of 6*3 irreducible representations are introduced; a 
doublet {Hi,H2)'^, a pseudo singlet Ha and a singlet Hs for each of up and down type Higgs 
sector. It is found from the tensor product (12. 2p that the most general supersymmetric Yukawa 
(mass) terms are written as 



where a, 6, ■ ■ ■ , / are independent Yukawa coupling constants. 

It is noticed that, if some of Higgs bosons have vanishing expectation values in the vacuum 
of the theory, corresponding mass matrix elements vanish and then mass (not Yukawa) texture 
is realized. That is the scheme we adopt in this paper for quark/lepton mass matrices. A non- 
trivial issue is whether such VEVs are obtained in the electroweak symmetry breaking vacuum. 
Therefore the Higgs potential should be carefully examined. 

Finally, the S'3-invariant bare Majorana mass for matter superfield \E'j {i = 1,2,3) is given 




Wy = ^Li {MD)ij ^R,- , 

aHi bHs + cHa dH2 



(2.3) 



(2.4) 



by 



m 



(2.5) 




(2.6) 
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3 5*3 Higgs scalar potential 

Let us consider supersymmetric theory with 5*3 flavor symmetry and introduce the foUowing 
Higgs superfields: 

Hus, Hds ■■ S3 singlets 

HuA, '■ S3 pseudo singlets . 

{Hy,i,Hu2), {Hdi,Hd2) ■ 5*3 doublets 

The symmetry-invariant supersymmetric Lagrangian is 



£ = y d^9d^9 K + l^j d^e W + h.c.j , (3.1) 

^ = E (^Le^"^na + ^ie^'^^da) , (3-2) 
a=S,A,l,2 

W — usHuS^Hds + IJ'AHuA^-HdA + AH) \ HuieHd,2 + Hu2^Hdi ) , (3.3) 



where Gu,d = ^^gyVy + (72^2 with the electrowcak vector superfields Vy and V2, and gy and 
g2 are the corresponding gauge coupling constants. The mass parameters ij>s,a,d are generally 
complex and e is the antisymmetric tensor for SU{2) weak indices (e^^ = 1). The Higgs scalar 
potential is then given by 

V = Kusy + Koft , (3.4) 

a=l,2,3 

+ \i,sf {hIsHuS + HlsHds) + {hI^H^a + hI^H^a) 

+ liinf (hIh^, + HlHar + hI^H^^ + Hl^Hd2) , (3.5) 

V^oft = '^IsHIsHuS + mlsHlsHds + {bsHy^seHas + h.c.) 

+'^I.aHIaHuA + mlj^Hl^HdA + {bAHy^AeHdA + h.c.) 
+mlj,(Hi,H^, + Hi,H^,) + mUHlHd, + H\^Ha2) 
+ [bD{H^ieHd2 + Hu2eHai)+h.c.] , (3.6) 
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where bx, mux and rudx {x — S, A, D) are the holomorphic and non-holomorphic mass parame- 
ters of supersymmetry breaking, respectively. The D terms are exphcitly given by 

Dy = --9Y [HlgHus - HlgHds + hI^H^a - H^^H^a 

+ hIiHui - H\^Hdi + HI2HU2 - H\2Hd2j , (3.7) 
= -g2{HlsT''Hus + H\sT^Has + Hl/r''HuA + H\j{r-HdA 

+ Hl^T'^Hui + Hl^T-Hd^ + Hl2T''Hu2 + i^la^'^^ds) , (3.8) 

with (a = 1, 2, 3) is the SU (2) weak generators. It should be noted that V^oft is introduced as 
the most general supersymmetry-breaking Higgs soft terms compatible with the 5*3 symmetry. 

We here mention a possibility of spontaneous breakdown of CP symmetry. The following 
terms in V^oft which contain complex couplings could be sources of CP violation: 

bsHuseHds + bAHuAeHdA + bo {HuieHd2 + Hu2eHdi) + h.c. . (3.9) 

Analyzing the stationary conditions for the Higgs potential, it is understood that the phases of 
Higgs VEVs should satisfy 9ua + ^da — ^ {a — S,A,1, 2) where 9ua{da) are the up(down)-type 
Higgs phases. While these phases do not directly appear in the Higgs scalar potential at the 
vacuum, they cannot be removed in general by field redefinitions unlike in the minimal super- 
symmetric standard model (MSSM). Therefore there is a possibihty of causing spontaneous CP 
breakdown by taking account of complex phases of Higgs VEVs in Higgs-mediated processes. 

4 Analysis of Higgs scalar potential 

Hereafter we assume that the charged Higgs fields do not develop nonzero VEVs so as to make 

the U{1)em symmetry remains intact after the electroweak symmetry breaking. The S3 Higgs 
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scalar potential we analyze thus becomes 

y = {\f^sf + mis) vls+ {If^sf + m%) vis - "^bsVusVds 

+ (I/zaI^ + ml^) vIj^ + (l/i^l V - 2hAVuAVdA 

+ + mlj^) [vl^ + vl^) + (l/iDl^ + mli^) [vl^ + vl^) 

-2hD {vuiVd2 + v^2Vdi) + ^^1±^X2, (4.1) 

o 

X = vis- + - + vli - vj^ + vl^ - vj^, (4.2) 
where Vx are the absolute values of Higgs scalars: 

V^S = \{H:s)\ , VuA = \{H:^)\ , Vul = KOI , Vu2 = \{H'^2)\ , 

Vds=\{Hl,)\, VdA = \{H%)\, Vdr = \{H',,)\ , Vd2 = \{H'd2)\ , (4-3) 

and the parameters hs^A,D have been chosen to be real positive by field redefinitions. Notice 
that the scalar potential (14. ip is invariant under the label-exchanging transformations, 1^2, 
S ^ A and/or u ^ d. 

First we discuss the instability of scalar potential at the origin of field space. The poten- 
tial parameters have to satisfy at least one of the following conditions in order to break the 
electroweak gauge symmetry: 

hi > {M' + mls){\fisf + mls) , (4.4) 
bl > (l/iAl' + m^Jd/UAl' + m^^) , (4.5) 
bl > i\f^Df + ml^){\fiDf + ml^). (4.6) 

Then the equations of motion become 

(lUsf + 'mls)vus = bsVds - X'vus , ill^sf + 'mls)vds = bsVus + X'vds , (4.7) 

{\^,A\' + mlj^)VuA = hAVdA - X'VuA , (IfiA^ + ml^)VdA = bAVuA + X'vdA , (4.8) 

+ ml^j)vui = bDVd2- X'vui , {\^,d\^ + mlj^)vd2 = bDVui + X'vd2 , (4.9) 
ilfj^nf + mljj)vu2 = boVdi - X'vu2 , (I/xdI^ + m^^)fdi = bDVu2 + A'vdi , (4.10) 



where X' = ^^-^X. These equations fl4.7p - fl4.lUp depend respectively on VuS{dS)i iJuA{dA)i '^Mi(d2) 
and fM2(di), except in the X' parts. The Z boson mass is given by 

^1 = \^9l + gl) {vis + vis + vIa + vIa + < + < + + vl^) . (4.11) 

In addition, substituting fl4.7p - fl4.10p for the scalar potential, we find the depth of the potential 
at the minimum: 

K.in = -^^^X^ < . (4.12) 
o 

This expression also holds in the case of 5*3 symmetry breaking model discussed in later section. 

Next we study the possibility of having vanishing VEVs in the minimum of the potential. 
In general, (14. 70 - 04. lOp are the coupled equations through the X parts. We separate these 
equations into three parts for the singlet (14. 7p . the pseudo singlet (14.80 . and the doublet (14.90 
and f l4.10p . As for the singlet and pseudo singlet parts, each coefficients of X' in (14.70 and (14.80 
are Vu{d)s aiid Vu(d)Ai respectively. Therefore vanishing VEVs makes the equations trivial within 
each sector. There is however a bit difference in the 5*3 doublet part [Eqs. (14. 9 p and (14.100 ], that 
is, they cannot be separated from each other. These equations have not only a relation through 
the X part but also a common parameter 6/) which originates from the symmetry invariance. 

4.1 5*3 singlet Higgses 

First let us examine the vanishing VEVs of 5*3 singlet Higgs fields. In this case, we take the 
other Higgs VEVs arbitrary. The analysis is the same for pseudo singlet Higgs fields. There 
are two possible patterns of zero VEVs; one is that both of the up and down-type Higgs VEVs 
are zero and another is that only one of them is zero. 

• VuS =VdS = ^ 

This solution always exists. Though such a solution does not make sense in models with 
one pair of Higgs doublets like the MSSM, it is possible in the present long as 

some of instability conditions for other Higgs parts are satisfied. 
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• Vus = 0, Vds 7^ or Vds = 0, Vus 7^ 

It is found that the parameter condition 65 = is necessary for this solution. For example, 
if one supposes Vus = and Vds 7^ 0, the stationary conditions mean 

2 2 

bs = and |/i5|' + mis ' ^-^^ {-^% + ^0) = 0, (4.13) 

where we have defined Xq = X\y^g=y^g=(). Therefore as long as the parameter bs is 
nonzero, the solution with one of VEVs being zero does not exist. It is noted that the 
condition 65 = is difficult to be satisfied exactly at any scale, and hence such a type of 
solution may not be realistic. 

4.2 5*3 doublet Higgses 

Including S3 doublet Higgs fields makes the theory complex but is promising. Unlike the case 
of ^3 singlet Higgses, all the four VEVs, Vui, Vu2, Vdi and Vd2 should be simultaneously taken 
into account because 6/5 is a common parameter as we mentioned above. In the analysis, we 
take the singlet Higgs VEVs, Vus, Vds, VuA and VdA, to be arbitrary. It is important to notice 
that the stationary conditions for S3 doublet Higgses are described (except in the X parts) by 
two pairs of VEVs: Eq. (14.91) for (f„i,fd2) and Eq. (14.101) for {vu2,Vdi). It can be seen from 
the result in the previous section that the existence of VEV pairs means that the VEV forms 
such as Vui = and Vd2 7^ are not allowed, i.e. only one VEV in each pair does not develop a 
nonzero VEV unless an unlikely condition bo = is satisfied. 

There are 16 (= 2^) patterns for the 4 VEVs of 53 doublets. Among them, only 7 patterns 
are theoretically independent due to the label-exchanging invariances of the potential: 1^2 
and u ^ d, mentioned before. For example, in the case that Vui = Vd2 = and Vu2, Vdi 7^ 0, 
the equations of motion give 

2\^^D\' + ml^ + = f ^ + (4.14) 

A similar form of vacuum equation is obtained for a theoretically equivalent case that Vu2 = 
Vdi = and Vui, Vd2 7^ by use of label exchanges. Therefore it is enough to consider either 
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Table 2: All representative VEV patterns of the 6*3 doublet Higgs fields. The blank entries 
denote no needs of parameter conditions. 
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Table 3: All VEV patterns of the S3 doublet Higgs fields without parameter conditions. 

of these patterns. All 7 possible patterns obtained from the potential analysis are shown 
in Table [2l It is noted that the parameter conditions given in Table [21 i.e. b£, = and 
2|/iDp + + ^dD ~ unlikely to be satisfied exactly, because these relations are not 

protected by symmetry arguments and generally unstable against quantum corrections (though 
there are some models for a vanishing b parameter at a given scale [IH])- Therefore in this paper 
we do not consider the solutions with non-trivial parameter conditions. In the end, we have 
found the relevant vacuum solutions without parameter conditions, which are presented in 
Table El 

To summarize the results in Section 4, we have found all possible minima of the scalar 
potential for 5*3 singlet and doublet Higgs fields (Table H]). These solutions do not require any 
exact tuning of Lagrangian parameters for electroweak symmetry breaking and are physically 
available. The potential depth at each minimum is controlled by Higgs mass parameters in the 
Lagrangian, and any of the VEV patterns in Table H] can be made the global vacuum of the 
theory, as seen in Eq. (14.12p . 
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Table 4: All possible minima of the scalar potential for 5*3 singlet and doublet Higgs fields 
without tuning of Lagrangian parameters for electroweak symmetry breaking. The blank entries 
denote non-vanishing VEVs. 

5 Quark and lepton mass textures 

As a phenomenological application of the vacuum analysis performed in the previous section, 
we study in this section the mass matrix forms of quarks and leptons derived from 5*3 flavor 
symmetry. Let us consider a supersymmetric theory with three-generation quark and lepton 
superfields in addition to the Higgs content previously analyzed. The three generations belong 
to a reducible 3 representation of S3, that is, two of three generations behave as a doublet and 
the other is a singlet or a pseudo singlet. When the first two generation superfields consist of 
5*3 doublets, the most general form of mass matrix is given by (12. 4p with non-vanishing Higgs 
VEVs. 

A usual approach often seen in the literature to obtain realistic forms of mass matrices 
is to control coupling constants of effective Yukawa operators with additional symmetries or 
to adjust coupling constants by hand so that the experimental data is reproduced. Unlike 
these approaches, our strategy in this paper is to dynamically realize mass matrix forms (mass 
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textures) in the vacuum of the theory. That is, as discussed in the previous section, some of 
Higgs VEVs vanish at the minimum of scalar potential, which in turn lead to texture (zero) 
forms of quark/lepton mass matrices. The available candidates of vacua are listed in Table HI 
One can see from the table that the VEV structures of up and down-type Higgs fields are not 
parallel and rather different due to the group properties of S3 flavor symmetry. This fact is 
favorable in light of the experimental data of fermion masses and mixing angles, which data is 
well known to show that the up and down quark sectors would have highly different generation 
structures. 

In the following, we assume as an example that the S3 representations of three-generation 
matter flelds are 2 -|- Is- If one adopts 1a instead of Is, phenomenological results are changed 
according to Eq. (12. 2p . Examining the vacuum patterns in Table HI we flnd that the following 
case with four zero VEVs leads to the simplest texture (i.e. the maximal number of zero matrix 
elements) with non-trivial flavor mixing: 



The VEV pattern obtained by interchanging the labels 1 ^ 2 is also the case. That can be 
covered by exhausting the generation label exchanges of matter flelds, and therefore we safely 
focus on the vacuum (15.11) and examine all types of 5*3 charge assignments of matter flelds. In 
what follows, we use the notation that Qi, Ui, di, Li, Cj, and z/j {i = 1, 2, 3) are the superflelds 
of left-handed quarks, right-handed up quarks, right-handed down quarks, left-handed leptons, 
right-handed charged leptons and right-handed neutrinos, respectively. 

5.1 Quark mass textures 

When the flrst and second generation quarks behave as 5*3 doublets, the up and down quark 
mass textures in the vacuum (15.11) are read from the generic form of mass matrix (12. 4p : 



Vul = Vd2 = VuS = VdS = 0, 



Vu2, Vdl, v.. 



,A, VdA 7^ 0. 



(5.1) 




(5.2) 
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where the blank entries mean zeros, and a, /5, ■ ■ ■ , 7 and 5 are the Yukawa couphng constants. 
In this flavor basis, the above matrices appear not to lead to appropriate mass hierarchies and 
generation mixing. We have exhausted the S3 charge assignments of matter fields and found 
that only the following 4 cases are almost consistent with the current experimental data of 
fermion masses and mixing angles. 

5.1.1 (Q2,Ql) + Q3, (U3,U2)+U1, (d2,d3)+di 

The title of this subsection means the S3 charge assignment that the second and first generations 
of left-handed quark superfields consist of a doublet and the third one Q3 is a singlet, and 
similarly for Ui and di. In this case, the quark mass matrices are given by 

M„ = I I , Md = \ ea bd \ . (5.3) 




In the first order approximation, they predict the quark mass eigenvalues and mixing angles 

rriu = , rric = Cu , rrit = iu , (5.4) 

ma = ad , ttIs = ea , mt = id , (5.5) 

= - ^ , V,b = ^-^ , Kfe = ^ - ^ , (5.6) 

which lead to a relation among the observables: 

- — = f + —V^ . (5.7) 

Such parameter-independent relation may be useful to examine whether the model can well 
describe the observations. If one evaluates the relation (15. 7p with respect to the mixing matrix 
element Vch, one third of the observed value is reproduced, which implies a necessity of some 
modification of the mass matrix forms. 
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5.1.2 (Q3,Q2) + Ql, (U1,U2)+U3, (d3,d2)+di 

In the second case, ((53,(52), (^1,^2) and (^3,^2) are the 5*3 doublets. The mass matrices then 
take the forms: 



= \ -du Cu , Md= \ dd -fd . (5 





The flavor charge assignments of Qi and di are the same, which leads to a generation-symmetric 
texture form like the Fritzsch ansatz [9J for the down quark mass matrix. Also such charge 
assignment might be relevant to the flipped SU{5) unified models [19]. The prediction for the 
quark mass eigenvalues and mixing angles is given by 



/ , n 5 ""C, — \/du~^^U 5 — iu , (5-9) 
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rud = , = ^ , m, = Id, (5.10) 

= ^ , = -ft , Vut= '-4^ , (5.11) 



Jd "d 



in the first order approximation. It is noted that the leading contributions to the mixing angles 
mainly come from the down quark sector, and then satisfy the following two relations among 
the observables: 



\Vub\ = \Vus\[ — ) , \Vcb\ = J—. (5.12) 

These do not so largely deviate from the experimental data, but predict a bit large (small) 
value for the mixing matrix element Vcb (for the mass eigenvalue m^). 

5.1.3 (Q2,Q3) + Ql, (U2,U3) + U1, (d2,di)+d3 

The third pattern is the case that {Q2, Q3), {u2, u^) and (^2, di) are the S3 doublets, which lead 
to the mass textures: 



3/2 



M^= \ du fu \ , Md= \ dd ed . (5.13) 
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This type of charge assignment is compatible with SU{5) grand unification [20]. The masses 
and mixing angles are determined in the first order as 

rriu = -^-^^ , mc = f , rrit = iu , (5.14) 



= Lf' 2 ' = + e\ , mfc = id , (5.15) 

Vus = , Vcb = , Vub= . (5.16) 

Jv. *u /" 

It is noted that the leading contributions to the mixing angles mainly come from the up quark 
sector, and then satisfy the following two relations among the observables: 

\Vcb\ = J— , \VusV,b\ = \Vub\ , (5.17) 



nit 

with which we would have some discrepancy between the prediction and the experimental data. 

5.1.4 (Q2,Q3) + Ql, (U2,U3) + U1, (di,d2)+d3 

The last case is the 5*3 charges; {Q2, Q3), (m2, U3) and {di, 1^2) are the doublets, which assignment 
is also consistent with SU{5) unification. Note that this type of assignment differs from the 
third case above, only for the doublet constitution of right-handed down quarks in the first and 
second generations. Therefore the third and this fourth cases are physically equivalent if all 
parameter spaces were taken into account. That is, these two cases correspond to two separate 
parameter regions in a single theory, both of which regions are phenomenologically viable. The 
mass texture forms now become 





Mu= \ du fu \ , Md= \ dd ed I , (5.18) 
\ -fu iu / \ -ed id ^ 

which give the masses and mixing angles in the first order approximation: 

rriu = -^^^ , = f , rrit = in , (5.19) 

md = , m, = Crf , rUb = id , (5.20) 

= -^ + ^, = -i^, = (5.21) 
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There are 8 free parameters for 9 observables, and we have one predictive relation 

\Vd,\ = . (5.22) 

V 

This relation is known to well fit the experimentally observed values and has already been 
discussed in other theoretical frameworks |16l |2T]. However the precise numerical estima- 
tion indicates that the equation f l5.22p is not exactly satisfied with the observational data: 
\Vch\ = 0.039 - 0.044 and \/m^lmt = 0.057 - 0.064 [22]. Some remedies can easily be found. 
The discrepancy is removed with radiative corrections, for instance, the renormalization-group 
effects on masses and mixing angles. A naive and probable source of such effects is the large 
top-quark Yukawa coupling. (The strong gauge coupling does not affect mass ratios and gen- 
eration mixing due to its flavor universality.) For example, in the MSSM we obtain the 
renormalization-group equation 



din [\Vcb\/ ^J^m^fmt 



^ [vl-vD- (5.23) 



(iln/i 327r2 

where [i is the renormalization scale and are the top and bottom Yukawa couplings, re- 
spectively. The positive coefficient of the top- Yukawa contribution implies that the discrepancy 
is indeed reduced in lower energy regime. It is however noted that in the present model the 
flavor symmetry is supposed to be broken at low energy such as the electroweak scale and hence 
the scale dependence might not be enough to make the relation fl5.22p fulfilled. 

Another possible source of radiative corrections comes from supersymmetry-breaking pa- 
rameters. The scalar fermions propagate in the loops and the chirality is fiipped via holomorphic 
or non-holomorphic couplings of scalar quarks. In both of these cases, the corrections do not 
change zero structure of mass textures but could modify nonzero matrix elements, depending 
on supersymmetry-breaking mass parameters of scalar quarks. For example, if the (M„)23,32 
elements receive such corrections, the relation fl5.22p is modified and the theory would become 
viable in light of the current experimental data. We leave a detailed study of model dependence 
on supersymmetry-breaking parameters of matter fields to future investigations. 
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In the present analysis, we have neglected CP- violating phases. All the mass matrix elements 
in this subsection can be taken to be real by quark phase redefinitions. This fact implies that 
the mass matrices are diagonalized by real orthogonal matrices up to overall phase rotations. 
In this case the observable CP phase is induced by the presence of different phase rotations 
between the up and down left-handed quarks. This is the case for the quark mass matrices in 
this subsection. The prediction of quark sector CP violation is consistent with the experimental 
data as long as the mixing angles are properly reproduced. 

5.2 Lepton mass textures 

Under the standard model gauge symmetry, the flavor charge assignment and Yukawa couplings 
of the lepton sector are free from those of the quark sector. One could explore the patterns 
of ^3 charges for lepton flelds, just as in the previous analysis for the quark sector, so that 
various phenomenological constraints are satisfled. However an attractive way to determine 
lepton flavor charges is to promote the theory to be embedded into grand uniflcation, where 
lepton and quark multiplets are unifled and have the same flavor charges. Along this line of 
thought, the cases 2, 3, and 4 in the previous section are the candidates to be considered. It is 
however easily found that the cases 2 and 3 are not reconciled with the experimental data of 
lepton sector, even if one includes additional Higgs flelds in higher- dimensional representations 
of unifled gauge symmetry. Thus we flnd the unique solution, the case 4 (in Section [5.1.41) . for 
5*3 flavor charges of lepton flelds which would be compatible both with uniflcation hypothesis 
and the observed data. It may be interesting to remind that, only from the analysis of quark 
sector performed in the previous sections, the case 4 has been phenomenologically singled out. 

Let us flrst consider the charged lepton sector. We assume here SU{5) grand uniflcation 
because the quark flavor charges in the case 4 are consistent with it. The left-handed charged 
leptons Li are combined with the right-handed down quarks into unifled gauge multiplets (anti 
quintuplets), and therefore (Li, L2) consists of an 5*3 doublet and L3 a singlet. Similarly, (e2, 63) 
transforms as an S3 doublet and Ci as a singlet. The charged lepton mass matrix Mg is given 
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by the transpose of that of down quarks Md in (15.181) and becomes at low energy 

/ dd 3ed \ 

hd -3ed 



Me =i r 



V 



(5.24) 



id j 



where we have included a group-theoretical factor '—3' [23j in front of the element e^. Such 
a facto. o.gi.a.es .o. a Yulcawa coupling to mghe.C.cnsional Higgs field |c.g. 45-plct of 

5'f/(5)] □ and is known to make the mass eigenvalues of charged leptons well fitted to the 
observed values if one takes account of renormalization-group effects on the down quark Yukawa 
couplings from the strong SU (3) gauge sector [that has been effectively included as the factor 
r in (I5.24p ]. The mass matrices Md (15.181) and (I5.24p are found to satisfy the relations 

3mp 



(5.25) 



md Srus 

which lead to a better explanation for the mass eigenvalues than that without the group- 
theoretical factor. We have found that all the other patterns to include factors —3 are not 
consistent with the observation. Therefore the charged lepton mass texture (15.241) is the unique 
simplest solution with our vacuum aligned scheme in SU (5) grand unification with S3 flavor 
symmetry. 

The mixing matrix which rotates the left-handed charged leptons to diagonalize Me is 



( 



V 



dd £d 

3ed id 
_3£d 1 



dd 3ea \ 
3ed id ' 



(5.26) 



in the first order approximation. It is found from this expression that |(t/e)i3| = "tTT-^i/^r — 1/17 
and |(t/e)23| = I(f^e)i2(^e)i3| ^ ^fi/^r- In particular, the latter means that the observed 
large lepton mixing between the second and third generations must come from the neutrino 
sector. As for (f/e)i2, a naive upper bound is obtained from the experimental upper bound on 
the 1-3 lepton mixing Ves; |(f/e)i2| < v^(K3)max5 where Vai (a = e,fi,T and i = 1,2,3) is 
the observable lepton mixing matrix [24J. Interestingly, a lower bound on the charged lepton 



^The representation of Higgs fields under unified gauge symmetry is independent of the potential analysis at 
the electroweak scale given in the previous section. 
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contribution to the 1-2 lepton mixing |(f/e)i2| can be expressed in terms of the observables by 
imagining that there should not be any fine tuning of parameters in reproducing V^s in ( I5.2ip . 
That gives |(t/e)i2| > 3(me/m^)/Ks ~ 1/16. 

We introduce the three generation right-handed neutrinos and utihze the seesaw mecha- 
nism [25] to obtain tiny neutrino masses. The neutrino Dirac mass texture in the vacuum 
of theory and the right-handed neutrino bare Majorana mass matrix M/j are respectively read 
from the fiavor symmetry invariance and given by 

/ K cA / Ml \ 

= l-b, , Mr= \M, . (5.27) 

\9. J \ M,J 

Here we have taken the charge of right-handed neutrinos as (z/i, z/2) + 1^3, that is, the first 
two generations make a doublet. It is however noticed that the charge assignment of z/j (equiv- 
alently the label changing effect of z/j) is completely irrelevant to low-energy physics and the 
generation structure of light neutrino mass matrix is determined only by the flavor charge of 
left-handed leptons Lj. After integrating out heavy modes, the light neutrino mass matrix 
Ml = -M^M^^M^ becomes 



R ^"1/ 

/ —z 1 —X \ 

1 2y —xy 



Ml = — 



(5.28) 



y —X —xy J 



K K K 

Taking into account that the charged lepton sector has small generation mixing, it is found 
that this form of neutrino mass matrix is suitable for large generation mixing with the inverted 
mass hierarchy of light neutrinos: m2 — mi ^ m^. The 2-3 large lepton mixing is controlled 
by X ~ 0{1). The observed mass squared differences imply the traceless condition 2y ~ z 
{\y\, \z\ ^ 1). In the limit y,z ^ 0, the 1-2 and 2-3 neutrino generation mixings are maximal 
{6'^2 — ^23 ~ '^/^) ^^"^ angle becomes zero {6'(^ = 0). Including flnite effects of y and z, 

we obtain the flrst order expressions: 

ia,Y? 9^n = 1 ,^ , tan^nq = x , tan^^'o = ^ , (5.30) 
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mi = ^(v/rT^-Z/ + |) , m2 = ^(v/rT^ + l/-|) , m3 = 0. (5.31) 



If the small mixing contributions from Me are neglected, there is a prediction which relates the 
observables in neutrino oscillation experiments: 

2tan^^5'o ^„ , ml — m?, 

+ tan - 1 = -r^o W • 5.32 

tan 4(m2-m2) ^ ^ 

This is certainly consistent with the current experimental data. Finally, including the charged 
lepton contribution f l5.26p . we obtain 



Ve2 


= sm9'(2 


-(f/e)i2COS^r2 


COS 023, + iUe)l3 COS 6^12 siu 6^23 , 


(5.33) 




= sin 6123 






(5.34) 




= sin ^5^3 


-(f/e)l2Sin^2'3 


- (f/e)l3 COS 6^23 , 


(5.35) 



where Oy^ <^ 1 has been used. It is found that Ves has a lower bound; \Ve3\ > 0.04, which is 
derived by taking account of the experimental data, Ve2, V^3, Am^2 ^'^is with the 
relation f l5.32p and the constraints on (f/e) 12,13 discussed above. That will be tested in future 
experiments such as the double Chooz [26]. 

In Section 5, we have shown that our scheme for generating zero texture forms with vacuum 
alignment is applied to mass matrices of quarks and leptons. Exhausting the patterns of matter 
flavor charges, we have found the highly limited numbers of charge assignments are phenomeno- 
logically viable. In particular, our scheme is also consistent with SU{5) grand unification, and 
the uniquely determined matter flavor charges predict typical low-energy phenomenology such 
as the relations among the observables, independently of model parameters. 



6 Higgs mass spectrum and 5*3 breaking terms 

The analysis in the previous section has shown that the vacuum (15. ip is phenomenologically 
interesting if applied to the quark and lepton mass textures. However the model has a problem 
that some massless bosons generally appear at any vacuum where more than one sets of Higgs 
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Hus 


HdS HdA 


Hui 


Hu2 


Hdi 


Hd2 


U{l)x 


+1 


-10 


+1 








-1 


U{l)y 


+ 1 


-1 +1 -1 


+ 1 


+1 


-1 


-1 



Table 5: The U{1) symmetries of the Higgs scalar potential with the 5*3 breaking terms (16. ip . 

bosons have non-vanishing VEVs. This general feature comes from the fact that the Higgs 
scalar potential given in Section 3 has an enhanced global symmetry SU (2) x f/(l)^ and leads to 
massless Nambu-Goldstone bosons in the electroweak broken phases. It is therefore reasonable 
to (softly) break the flavor symmetry within the scalar potential so as not to enhance the 
global symmetry. In this paper we introduce supersymmetry-breaking soft terms which violate 
5*3 flavor symmetry. We flnd that the following simple set of breaking terms is appropriate for 
the present purpose: 

= bsDHusHd2 + KsoHuiHds + bAoHuAHdi + b'j^j;)Hu2HdA + ^-c. , (6.1) 

There are several reasons why these terms are chosen: (i) The global symmetry of total Higgs 
potential is broken down to U{l)x under which H^s and Hui {Hds and Hd2) are positively 
(negatively) charged. Note that, combined with the f/(l)y hypercharge symmetry, U{l)x is 
converted to another global U{1) where Hua and Hdi {HdA and Hu2) are positively (negatively) 
charged (see Table [5]). As discussed in Section |42| there are four sets of up and down type Higgs 
fields in the original potential without flavor-breaking terms, which potential therefore has four 
f/(l) invariances. Each set of {Hu,Hd) are vector-like fields under these f/(l)s. Including the 
soft terms (16.11) economically breaks two of them, and hence U{1)x,y in Table [5] remain intact. 
It is important to notice that, in the electroweak broken vacuum (15. ip . U{l)x is unbroken 
and the Nambu-Goldstone boson of f/(l)y is partly absorbed by the massive gauge boson Z 
and becomes unphysical. Therefore no physical Nambu-Goldstone fields are associated with 
symmetry breaking. We numerically checked by examining the neutral Higgs mass matrix that 
there appear no massless scalar bosons, (ii) The second reason is that, when including the 
breaking terms (16. ip . the stationary conditions can still be solved without imposing any tuning 
of model parameters, (iii) Flavor breaking effects in supersymmetry-breaking holomorphic mass 
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terms do not propagate to other sectors. Such a favorable property for the model simplicity 
does not hold in cases that flavor symmetry violation resides in other supersymmetry-breaking 
sectors. 

The model also predicts a light Higgs boson exactly parallel to the MSSM case. In order to 
find an approximate form of the Higgs mass spectrum, we take a simplifying assumption: 



I |2 I 2 



I |2 I 2 



— 2 



[x = S, A, D) 



JSD 



b'sD = ^AD 



AD 



(6.2) 

6' , (6.3) 

with a hierarchy m^, \b'\ ^ f ^ = t>^2 + v'^i + "^uA + '^dA- this case, we can analytically 
write down the mass matrices of Higgs fields. In particular, for the real parts of neutral Higgs 
bosons, the mass matrix is given by 



Mi 



"'ul 

Ka 

Ks 
K2 

'^dA 
'''dS 



( vr? 



l)2di _|_ Jj'^MA 

Vu2 Vu2 



-b 



jjMA. _|_ y ""di 



-b 



-b 



rn? —b' 
—V fn? 



-b 
-b' 



^ Vu2 _j_ yVuA 



-b' 
-b 



-b' \ 



-b 



\ 



+ b' 



I Vu2 
VdA 



— 2 

m 



J 

(6.4) 



where {x = m1,m2, ■■■) are the neutral components of Higgs bosons. Diagonalizing this 
matrix, we obtain the neutral Higgs masses squared: 

{M^o, M^o, Mlo, M^o, Mlo} = {0{v^), m^-b-b', m^+b-b', m^-b+b', m^+b+b'}, (6.5) 

and the other three mass eigenvalues squared are written in rather complicated expressions, 
but of the order of 0{m'^,b,b'). We therefore have a light mode with a weak scale mass 0{v) 
and its eigenfunction is explicitly given by 



= - {VuAhlA + VdAh% + ^^«2/i°2 + Vdlhdi) ■ 
V 



(6.6) 



This light Higgs mode receives a sizable radiative correction from large top Yukawa coupling 
similar to the MSSM, and would be made heavy enough to satisfy the lower bound from the 
LEP experiment [22]. 
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Figure 1: A Higgs- mediated tree- level FCNC process for the K meson system. 

7 Tree-level FCNC 

Since there are multiple electroweak doublet Higgses which couple to matter fields, fiavor- 
changing processes are mediated at classical level by these Higgs fields. In the previous section, 
we show that all but one Higgs bosons have masses of the order of supersymmetry-breaking 
parameters. Therefore the experimental (un)observations of FCNC rare events would lead to a 
bound on the supersymmetry breaking scale. Among various experimental constraints, we find 
the most important constraint comes from the neutral K meson mixing, which gives a lower 
bound on heavy Higgs masses being larger than a few TeV. 

The down quark mass matrix in Eq. ( 15.18^ is diagonalized as 



UlMdU^ 



rrib 



Ur 



/ 1 

v 



bddd eddd 



hi 

eddd 



\ 



J 



u, 



R 



1 


dd. 






ed 


id 


_dd 


1 





ed 




ed 


dd 


1 


id 


id 





(7.1) 



(7.2) 



Here and hereafter in this section, we simply take ~ d^. The other parameters are fixed by 
the mass eigenvalues. Tree-level FCNC is mediated by Higgs fields like the diagram shown in 
Fig. [H For the heavy mass eigenstates [i = 1, 2, 3, 4), the tree-level K^-Ks mass difference 
^^tree jg gjygj^ ^]-^g matrix element of the effective Hamiltonian between K mesons j27] . 
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which is analytically evaluated from the above mixing matrices in the present 



\2/ \2 ( / \ 2 T ^ ' 

\ / m \ _l/ m \ II 

, (7.3) 



K 



rriK \^ ( ms\^ "sfXi V I ^ 



ms + mdj \mb J \\ms + mdj 2 

where rriK and fx are the mass and decay constant of K meson and Mh is an average of the 
Higgs masses, l/Mfj = l{l/Mlo + l/M^o + + l/M^o). The parameter r/ contains the 

down-type quark Yukawa couplings which do not contribute to mass terms and explicitly given 

by 

^ _ iyd)22bdVdi {yDisddVdi 

ml msmb ' 

where denotes the matrix of Yukawa coupling of down-type quarks to the Ss singlet Higgs 
boson. The 5*3 singlet Higgs Hs has a vanishing VEV in the present vacuum, and therefore r] 
is regarded as a free parameter. For example, r/ = if one does not include Hs in the theory. 
The other Higgs fields also contribute to the K meson mixing in a similar order. 

In order to estimate the bound on supersymmetry breaking scale, we calculate the ratio 
of the exact numerical value Am^°° to the standard contribution Am^. In Fig. [21 we show 
the ratio versus the averaged Higgs mass Mh for the cases of ?7 = and r] = 0.03 as typical 
values. It is found from the figure that heavy Higgs masses, in turn the supersymmetry breaking 
masses, are bounded from below so as to suppress the extra Higgs contribution compared with 
the standard model one, which bound is roughly given by 

f 3.8 TeV (ri = 0) 
Mh > < ' ' , (7.5) 

[ 1.4 TeV (r/ = 0.03) 

where we have used the experimental data rriK = 490 MeV, fx = 160 MeV and taken Vdi = 
100 GeV as a typical electroweak scale. Other flavor- changing rare processes are also suppressed 
for such heavy Higgs fields of few TeVs. For example, the /i — e -|- 7 process is given in terms 
of extra Higgs masses and lepton Yukawa couplings [28]. The branching ratio is found to be 
suppressed enough under the constraint of eq. (17.51) . 
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Figure 2: Am^^^'^/Am^^ as the function of Higgs mass parameter Mh- The sohd and dashed 
hues correspond to = and 0.03, respectively. 



8 Summary 

In this paper we have discussed the structure of Higgs potential and fermion mass matrices in 
supersymmetric models with ^3 flavor symmetry. The electroweak doublet Higgs fields belong 
to non-trivial representations of the flavor symmetry and lead to restricted forms of vacuum 
structure. We have examined possible zero elements (textures) of quark and lepton mass 
matrices. Our approach to have mass textures is rather different from previous ones in a sense 
that flavor symmetry does not forbid coupling constants of (effective) Yukawa operators, but 
leads to mass texture zeros by ensuring some of Higgs fields have vanishing expectation values 
in the vacuum of the theory. The vanishing mass matrix elements are obtained dynamically in 
the vacuum without any exact tuning of model parameters, and their positions are controlled 
by flavor symmetry. An interesting point is that, due to the flavor group structure, the up 
and down quark mass matrices are automatically made different, which lead to non- vanishing 
generation mixing. We have exhausted the patterns of flavor symmetry charges of matter fields 
and found that the simplest, viable forms of mass matrices including the neutrino sector are 
uniquely determined. In particular the lepton mixing Vgs is predicted within the range that 
will be tested in near future experiments. 

We have also discussed the physical mass spectrum of Higgs bosons and its phenomeno- 
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logical implication. The flavor symmetry is softly broken by a certain class of holomorphic 
supersymmetry-breaking mass terms of Higgs bosons. That neither destabilizes the desired 
vacuum nor introduces any fine tuning of coupling constants. Given these flavor breaking 
terms, it is found that all Higgs bosons, except for the lightest one, can be made heavier than 
a few TeV enough to satisfy the experimental bounds such as that from the neutral K meson 
mixing. Detail phenomenological analyses including CP violation and FCNC will be presented 
elsewhere. 

We would like to thank J. Kubo, T. Kobayashi and H. Nakano for helpful discussions. This 
work was supported in part by scientific grants from the Ministry of Education, Science, Sports, 
and Culture of Japan (No. 17540243 and 17740150), and grant-in-aid for the scientific research 
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discoveries of Higgs boson and supersymmetry" (No. 16081209). 



26 



References 

[1] S. Pakvasa and H. Sugawara, Phys. Lett. B 73 (1978) 61; 
D. Wyler, Phys. Rev. D 19 (1979) 330; 

L.J. Hall and H. Murayama, Phys. Rev. Lett. 75 (1995) 3985; 
R. Dermisek and S. Raby, Phys. Rev. D 62 (2000) 015007; 

R.N. Mohapatra, A. Perez-Lorenzana and C.A. de Sousa Pires, Phys. Lett. B 474 (2000) 

355; 

J. Kubo, A. Mondragon, M. Mondragon and E. Rodriguez- Jauregui, Prog. Theor. Phys. 

109 (2003) 795 [Erratum-ibid. 114 (2005) 287]; 

T. Kobayashi, J. Kubo and H. Terao, Phys. Lett. B 568 (2003) 83; 

P.P. Harrison and W.G. Scott, Phys. Lett. B 557 (2003) 76; 

J. Kubo, H. Okada and F. Sakamaki, Phys. Rev. D 70 (2004) 036007; 

S.-L. Chen, M. Frigcrio and E. Ma, Phys. Rev. D 70 (2004) 073008 [Erratum-ibid. D 70 

(2004) 079905]; 

W. Grimus and L. Lavoura, JHEP 0508 (2005) 013; 
Y. Koide, Phys. Rev. D 73 (2006) 057901. 

[2] N. Haba and K. Yoshioka, Nucl. Phys. B 739 (2006) 254. 

[3] Y. Yamanaka, H. Sugawara and S. Pakvasa, Phys. Rev. D 25 (1982) 1895 [Erratum-ibid. 

D 29 (1984) 2135]; 

C. Hagedorn, M. Lindner and R. N. Mohapatra, JHEP 0606 (2006) 042. 

[4] W. Grimus and L. Lavoura, Phys. Lett. B 572 (2003) 189; 

W. Grimus, A. S. Joshipura, S. Kaneko, L. Lavoura and M. Tanimoto, JHEP 0407 (2004) 
078; 

W. Grimus, A. S. Joshipura, S. Kaneko, L. Lavoura, H. Sawanaka and M. Tanimoto, Nucl. 
Phys. B 713 (2005) 151. 



27 



[5] K.S. Babu and J. Kubo, Phys. Rev. D 71 (2005) 056006; 

Y. Kajiyama, E. Itou and J. Kubo, Nucl. Phys. B 743 (2006) 74. 

[6] M. Frigerio, S. Kaneko, E. Ma and M. Tanimoto, Phys. Rev. D 71 (2005) 011901. 

[7] E. Ma and G. Rajasekaran, Phys. Rev. D 64 (2001) 113012; 
E. Ma, Mod. Phys. Lett. A 17 (2002) 2361; 

K.S. Babu, E. Ma and J.W.F. Valle, Phys. Lett. B 552 (2003) 207; 

M. Hirsch, J.C Romao, S. Skadhauge, J.W.F. Valle and A. Villanova del Moral, Phys. 

Rev. D 69 (2004) 093006; 

A. Zee, Phys. Lett. B 630 (2005) 58; 

G. Altarelh and F. Feruglio, Nucl. Phys. B 720 (2005) 64; 

X.-G. He, Y.-Y. Keum and R.R. Volkas, JHEP 0604 (2006) 039; 

E. Ma, H. Sawanaka and M. Tanimoto, jhep-ph/0606103 

[8] E. Derman and H.S. Tsao, Phys. Rev. D 20 (1979) 1207; 

D. Chang, W.Y. Keung and G. Senjanovic, Phys. Rev. D 42 (1990) 1599; 

D. B. Kaplan and M. Schmaltz, Phys. Rev. D 49 (1994) 3741; 

P.H. Frampton and T.W. Kephart, Int. J. Mod. Phys. A 10 (1995) 4689; 
A. Aranda, C.D. Carone and R.F. Lebed, Phys. Lett. B 474 ( 2000) 170; 
N. Haba, A. Watanabe and K. Yoshioka, Phys. Rev. Lett. 97 (2006) 041601; 
C. Hagedorn, M. Lindner and F. Plentinger, Phys. Rev. D 74 (2006) 025007. 

[9] H. Fritzsch, Phys. Lett. B 73 (1978) 317; 

P. Ramond, R.G. Roberts and G.G. Ross, Nucl. Phys. B 406 (1993) 19. 

[10] H. Nishiura, K. Matsuda and T. Fukuyama, Phys. Rev. D 60 (1999) 013006; 

E. K. Akhmedov, G.G. Branco and M.N. Rebelo, Phys. Rev. Lett. 84 (2000) 3535; 
S.K. Kang and G.S. Kim, Phys. Rev. D 63 (2001) 113010. 



28 



[11] P.H. Frampton, S.L. Glashow and D. Marfatia, Phys. Lett. B 536 (2002) 79; 
Z.-z. Xing, Phys. Lett. B 530 (2002) 159; 

R. Barbieri, T. Hambye and A. Romanino, JHEP 0303 (2003) 017; 
A. Ibarra and G.G. Ross, Phys. Lett. B 575 (2003) 279. 

[12] A. Kageyama, S. Kaneko, N. Shimoyama and M. Tanimoto, Phys. Lett. B 538 (2002) 96; 
S. Kaneko and M. Tanimoto, Phys. Lett. B 551 (2003) 127; 

M. Honda, S. Kaneko and M. Tanimoto, JHEP 0309 (2003) 028; Phys. Lett. B 593 (2004) 
165; 

S. Kaneko, H. Sawanaka and M. Tanimoto, JHEP 0508 (2005) 073. 

[13] M.-C. Chen and K.T. Mahanthappa, Phys. Rev. D 68 (2003) 017301; 
M. Bando and M. Obara, Prog. Theor. Phys. 109 (2003) 995; 
M. Bando, S. Kaneko, M. Obara and M. Tanimoto, Phys. Lett. B 580 (2004) 229. 

[14] L. Lavoura, Phys. Lett. B 609 (2005) 317; 

W. Grimus, S. Kaneko, L. Lavoura, H. Sawanaka and M. Tanimoto, JHEP 0601 (2006) 
110. 

[15] M. Frigerio and A. Y. Smirnov, Phys. Rev. D 67 (2003) 013007. 

[16] N. Uekusa, A. Watanabe and K. Yoshioka, Phys. Rev. D 71 (2005) 094024; 
A. Watanabe and K. Yoshioka, JHEP 0605 (2006) 044. 

[17] C. Hagedorn and W. Rodejohann, JHEP 0507 (2005) 034; 

A. Merle and W. Rodejohann, Phys. Rev. D 73 (2006) 073012. 

[18] M. Dine, Y. Nir and Y. Shirman, Phys. Rev. D 55 (1997) 1501; 
R. Rattazzi an U. Sarid, NucL Phys. B 501 (1997) 297; 
M. Yamaguchi and K. Yoshioka, Phys. Lett. B 543 (2002) 189. 

[19] J.P. Derendinger, J.E. Kim and D.V. Nanopoulos, Phys. Lett. B 139 (1984) 170; 

L Antoniadis, J.R. EUis, J.S. Hagehn and D.V. Nanopoulos, Phys. Lett. B 194 (1987) 231. 

29 



[20] H. Georgi and S.L. Glashow, Phys. Rev. Lett. 32 (1974) 438. 

[21] K.S. Babu and Q. Shafi, Phys. Rev. D 47 (1993) 11; 

K.S. Babu and R.N. Mohapatra, Phys. Rev. Lett. 74 (1995) 2418. 

[22] S. Eidelman et al. [Particle Data Group Collaboration], Phys. Lett. B 592 (2004) 1. 

[23] H. Georgi and C. Jarlskog, Phys. Lett. B 86 (1979) 297; 

S. Dimopoulos, L.J. Hall and S. Raby, Phys. Rev. Lett. 68 (1992) 1984. 

[24] Z. Maki, M. Nakagawa and S. Sakata, Prog. Theor. Phys. 28 (1962) 870. 

[25] T. Yanagida, in Proceedings of the "Workshop on the Unified Theory and the Baryon 
Number in the Universe", Tsukuba, Japan, 1979, edited by O. Sawada and A. Sugamoto, 
KEK Report No. KEK-79-18, p. 95; Prog. Theor. Phys. 64 (1980) 1103; 
P. Ramond, in Talk given at the Sanibel Symposium, Palm Coast, Fla., 1979, preprint 
CALT-68-709; 

S.L. Glashow, in Quarks and leptons, proceedings of the advanced study institute (Cargese, 

Corsica, 1979), J.-L. Basdevant et al. eds.. Plenum, New York 1981; 

for an early work, see P. Minkowski, Phys. Lett. B 67 (1977) 421; 

See also R.N. Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44 (1980) 912. 

[26] K. Anderson et al, "White paper report on using nuclear reactors to search for a value of 



^is" , hep-ex/0402041 



[27] F. Gabbiani, E. Gabrielh, A. Masiero and L. Silvestrini, Nucl. Phys. B 477 (1996) 321. 
[28] For example, E. O. Dtaqn, JHEP 0402 (2004) 065. 



30 



